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I. 


INTRODUCTION 


Axial  compressor  design  systems,  involving  iterations 
between  axi-symmetric  and  blade-to-blade  calculations  and  in¬ 
corporating  empirical  representations  for  losses,  have  been 
remarkably  successful  in  producing  high  performance  machines, 
even  when  the  internal  flow  conditions  are  transonic.  Exam¬ 
ples  exist  today  of  commercially  successful  compressors  with 
multiple  transonic  stages.  However,  the  situation  remains 
that  the  full  range  of  performance  of  an  entirely  new  compres¬ 
sor  of  advanced  design  can  not  be  predicted  reliably  before 
the  machine  is  tested,  nor  can  the  steady  and  unsteady  loads 
on  the  blades  be  predicted  with  sufficient  accuracy  to  guar¬ 
antee  flutter-free  operation  and  to  allow  the  use  of  minimum 
structural  weight.  As  a  consequence,  compressor  development 
follows  an  evolutionary  procedure  in  which  problems  are  dis¬ 
covered  and  remedied  through  testing. 

Clearly, analysis  programs  are  needed  which  adequately 
describe  the  detailed  behavior  of  the  flow  in  high  speed  com¬ 
pressor  stages.  The  difficulties  however  are  enormous;  the 
flow  field  involves  embedded  shock  waves  in  the  presence  of 
boundary  layers,  leakage  flows  and  adverse  pressure  gradients 
and  passes  from  rotating  to  stationary  blade  rows.  Since 
severe  th  «*  elemental  processes  are  still  lacking  a 
fundamental  description,  and  since  the  three  dimensional  and 
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unsteady  nature  of  the  flew  is  inherent,  it  is  quite  .mlikely 

that  an  analysis  which  is  totally  devoid  of  "modelling"  ap¬ 
proaches  will  be  possible  in  the  foreseeable  future.  Conse¬ 
quently,  in  order  to  develop  computational  codes  for  compressor 
analysis  which  are  both  realistic  and  useful,  it  is  of  para¬ 
mount  importance  to  determine  experimentally  and  concurrently 
the  actual  flow  field  in  stage  geometries  to  which  the  develop¬ 
ing  analyses  can  be  applied.  While  such  an  activity  has  some¬ 
times  been  labelled  "code  verification",  it  is  more  likely  in 
fact  that  some  experimental  information  will  be  needed  during 
the  process  of  code  development.  For  example,  the  mechanism 
responsible  for  the  unreasonably  high  losses  observed  in  the 
outer  30%  of  several  transonic  rotors  ^  is  unlikely  to  be 
resolved  by  purely  analytical  approaches. 

The  present  work  follows  the  above  arguments.  It  re¬ 
ports  the  first  results  in  a  program  aimed  at  the  determina¬ 
tion  of  the  flow  field  produced  by  a  single  compressor  rotor 

in  a  rotor-first  arrangement.  Computation  of  the  rotor  flow 

(2) 

using  recently  developed  codes  is  intended,  however  only 
the  experimental  program  is  reported  herein. 

A  small,  450  HP  transonic  single-stage  axial  compressor 
and  test  rig  designed  for  continuous  operation  has  been  used. 
The  rig  was  built  to  facilitate  detailed  measurements  of  in¬ 
ternal  flow  behavior  and  specifically  to  provide  the  means  to 
obtain  information  with  which  to  verify  and  help  improve 
methods  of  flow  corn'  atation.  The  present  stage  is  one  of 
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in-house  design. ^  Instrumentation  was  developed  and  applied 

first  to  measure  time-averaged  flow  behavior  through  the  com- 
( 4 ) 

pressor.  Then  concentration  was  placed  on  the  determination 

of  the  detailed  internal  aerodynamics.  Clearly,  in  order  to 
measure  the  unsteady  flow  field  one  must  either  use  high  re¬ 
sponse  probes  or  an  optical  technique  of  some  kind.  While 
non-intrusive  laser  velocimetry  techniques  had  many  advantages, 
it  was  clear  that  they  required  an  expert  to  apply  them  in 
each  new  situation,  and  they  were  normally  able  to  measure  only 
two  components  of  the  velocity.  They  did  not  provide  measure¬ 
ments  of  the  pressure  field,  and  in  highly  unsteady  flows  the 
accuracy  was  questioned.  There  was  good  reason  therefore  to 
approach  the  problem  of  measuring  the  unsteady  flow  field  in 
other  ways,  and  ideally  to  obtain  independent,  redundant  mea¬ 
surements  of  the  same  phenomena.  The  work  reported  here  was 
such  an  attempt. 

Synchronized  (phase-locked)  digital  sampling  from  Kulite 
transducers  was  adopted  as  a  technique  for  defining  both  case 
wall  pressure  signatures  and  the  exit  flow  field  from  the 
rotor.  Early  results  were  reported  and  a  method  was  proposed 
in  Ref.  5  for  obtaining  the  full  periodic  velocity  field  at 
the  rotor  exit  from  measurements  with  two  very  simple  Kulite 
probes.  The  present  paper  reports  the  first  distributions 
measured  of  the  velocity  vector  (magnitude,  yaw  angle  and  pitch 
angle)  across  two  rotor  blade  passages  in  both  hub-to-shroud 
and  blade-to-blade  directions,  made  using  the  proposed 
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gives  an  account  of  the  experimental  procedures  used.  Results 
are  presented  and  are  discussed  first  on  the  basis  of  their 
qualitative  behaviour,  and  then  for  their  quantitative  accuracy 
Since  measurement  redundancy  was  inherent  in  the  approach 
it  could  be  shown  definitively  that  the  measurements  obtained 
in  regions  of  the  flow  which  were  periodic  in  the  machine 
frame  (steady  with  respect  to  the  rotor)  were  reliable,  whereas 
measurements  obtained  in  the  unsteady  (wake)  regions  could  not 
be  interpreted  accurately  on  the  basis  of  probe  calibrations 
obtained  in  a  steady  flow.  The  full  implications  of  these 
observations  are  discussed  in  detail. 


II.  DPDS  MEASUREMENT  APPROACH 


Concept 

Two  Kulite  semi-conductor  pressure  probes  of  very  simple 
design  are  positioned  behind  the  compressor  rotor  as  3hown  in 
Fig.  1.  The  "Type  A"  and  "Type  B"  probes,  shown  in  Fig.  2,  are 
mounted  eccentrically  in  plugs  so  that  they  can  be  rotated 
about  the  centers  of  their  tips.  Data  from  each  probe  is  di¬ 
gitized  under  "Pacer"  control.  Inputs  to  the  Pacer  are  one- 
per -blade  and  one-per-rev.  signals  obtained  as  shown  in  Fig. 

3.  The  complete  data  system,  based  on  a  Hewlett-Packard  HP21MX 
computer,  is  shown  in  Fig.  4.  The  Pacer  allows  the  data  to  be 
converted  at  any  point  in  the  rotor's  rotation  by  internally 
generating  a  pulse  train  with  a  frequency  which  is  phase-locked 
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to  the  blade  passing  frequency  but  multiplied  by  a  factor  of 
128.  Conversion  is  controlled  to  occur  when  a  programmed 
number  of  counts  matches  the  number  of  pulses  counted  by  the 
Pacer  from  the  one-per-rev.  reset  pulse. 

The  Pacer  control  allows  data  to  be  digitized  from  each 
probe  when  it  is  at  precisely  the  same  location  with  respect 
to  the  rotor  blading.  In  effect,  two  differently  oriented 
sensors  (Type  A  and  Type  B)  can  be  sampled  at  the  same  point 
in  the  rotor  frame  of  reference.  However,  since  each  probe 
can  also  be  rotated  about  its  tip,  data  can  also  be  obtained 
(at  each  rotor  point,  blade-to-blade)  as  a  function  of  the 
sensor  yaw  angle.  In  practice,  data  are  taken  at  each  of  256 
points  across  adjacent  pairs  of  rotor  blade  passages  with  both 
probes  rotated  in  unison  to  from  5  to  9  separate  yaw  angle 
settings.  Multiple  samples  (from  10  to  40)  on  successive 
revolutions  are  taken  at  each  point  and  probe  angle  setting 
and  ensemble-averaged  before  recording.  An  example  of  the 
ensemble-averaged  data  taken  on  the  compressor  annulus  center- 
line  across  two  particular  blade  passages  with  a  Type  A  probe 
is  shown  in  Fig.  5.  A  composite  picture  of  the  complete  data 
obtained  from  the  Type  A  and  Type  B  probes  for  one  radial  lo¬ 
cation  across  two  rotor  blade  passages  is  shown  in  Fig.  6. 

The  goal  then  is  to  reduce  the  information  obtained  in 
Fig.  6  to  obtain  the  distribution  of  the  velocity  vector  in 
the  blade-to-blade  direction. 
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Data  Reduction  Methods 


The  reduction  to  velocity  requires  the  following 
assumptions: 

i)  that  the  probes  respond  with  negligible  error  to  fluc¬ 
tuations  in  the  flow  which  are  at  blade-passing 
frequency; 

ii)  that  the  flow  is,  on  the  average,  steady  in  the  rotor 
frame. 

The  first  assumption  is  required  since  the  calibration 
of  the  probes  can  only  be  carried  out  in  a  steady  flow  field. 
The  second  assumption  is  required  since  the  two  sensors  are  in 
different  peripheral  locations.  (The  validity  of  the  assump¬ 
tions  will  be  discussed  after  an  examination  of  the  results.) 

In  general,  at  each  point  in  the  blade-to-blade  direction 
there  are  five  unknown  quantities;  namely,  flow  yaw  angle,  flow 
pitch  angle,  flow  velocity,  pressure  and  temperature.  In  all 
approaches  followed  to  date,  the  probe's  pressure  response  has 
been  expressed  as  a  function  of  the  Mach  number  and  it  is  as¬ 
sumed  the  effect  of  Reynolds  number  between  the  probe  calibra¬ 
tion  and  application  is  small.  The  unknowns  at  each  point  are 
then  yaw  angle,  pitch  angle,  Mach  number  and  pressure.  In  re¬ 
ducing  Mach  number  to  velocity,  the  stagnation  temperature  is 
assumed  to  be  constant  across  the  blade  passage  at  the  time 
averaged  value  measured  at  the  rotor  exit.  Thus,  in  principle, 
only  four  independent  pressure  measurements  are  required  in  or¬ 
der  to  derive  the  properties  of  the  flow  at  each  point.  Using 
5  to  9  probe  angle  settings  with  two  probes,  10  to  18  indepen¬ 
dent  measurements  result,  so  that  considerable  redundancy  is 


13 


always  present.  lit  should  be  noted  that  the  experimental 
procedure  described  below,  includes  on-line  calibration  of 
the  transducers  to  remove  the  effect  of  their  temperature 
sensitivity. ) 

Three  fundamentally  different  approaches  to  the  data 
reduction  have  been  examined.  The  methods  are  best  under¬ 
stood  with  reference  to  Fig.  7  which  shows  an  example  of  data 
obtained  at  a  single  point  in  the  blade-to-blade  direction 
(for  example,  any  two  corresponding  curves  in  the  two  data 
sets  in  Fig.  6) .  It  is  assumed  for  the  moment  that  the  data 
are  known  as  pressures  after  the  transducer's  have  been  cali¬ 
brated,  on-line. 

Method  1.  As  reported  in  Ref.  6,  in  this  method  the 
calibration  curves  for  the  two  probes  are  approximated  by 
fourth  order  polynomials  and  hence  the  flow  yaw  angle  (indi¬ 
cated  in  Fig.  7) ,  and  values  of  pA  fflax  and  pB  max  can  be  found 
from  the  test  data  by  curve  fitting.  The  Mach  number,  pitch 
angle  and  static  pressure  are  then  calculated  using  three  ap¬ 
proximate  equations  to  which  the  calibration  data  have  pre¬ 
viously  been  reduced.  The  values  of  pA  maix ,  pQ  m3>x ,  pA  and  aA 
are  sufficient  to  determine  the  unknowns,  and  hence  each  data 
point  obtained  with  the  Type  A  probe  which  is  well  away  from 
the  maximum  point  determines  a  result  for  the  flow  vector. 

Method  2.  As  reported  in  Ref.  7,  this  is  a  numerical 
iteration  approach  in  which  the  calibration  data  are  stored 
in  arrays  and  any  three  (3)  values  from  the  Type  A  probe  and 
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one  (1)  from  the  Type  B  probe  can  be  used  to  derive  the  four 
unknowns.  In  principle  the  approach  can  be  used  if  the  out¬ 
put  of  the  probe  is  not  symmetric  about  its  maximum  so  long 
as  only  one  maximum  exists. 

Method  3.  Reported  in  Ref.  8,  this  method  is  similar  to 

that  originally  proposed  in  Ref.  S  but  makes  use  of  analytical 

C9) 

tools  developed  more  recently.  It  uses  the  values  pA  max# 

pQ  _ _  and  p,  taken  at  a  discrete  value  of  a,  (here  63°) ,  as 

if  they  were  pressures  from  a  conventional  multiple-sensor  flow 
probe.  (The  values  must  be  obtained  from  the  data  by  interpo¬ 
lation.)  It  is  shown  in  Ref.  10  that  the  calibration  of  multi¬ 
ple  sensor  probes  can  be  well  represented  as  two  surfaces,  one 
for  Mach  number,  one  for  pitch  angle  as  a  function  of  two  pres¬ 
sure  coefficients  defined  in  terms  of  the  four  individual  sen¬ 
sor  measurements.  In  the  present  application,  the  flow  yaw 
angle  is  first  derived  by  a  curve  fit  to  the  Type  B  probe  data 
and  the  value  of  pA  is  taken  where  the  right  and  left  hand  sides 
of  the  Type  A  characteristic  are  separated  by  126  degrees. 

The  results  presented  in  the  present  paper  were  reduced 
using  Method  3  only. 


III.  EXPERIMENTAL  PROGRAM 
Summary 

The  measurements  reported  here  were  made  with  two  gen¬ 
erations  of  probes.  First,  probes  having  the  tip  geometry 
shown  in  Fig.  2  were  calibrated  in  a  free- jet  and  applied  in 
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the  compressor.  The  probe  tips  were  2.36  mm  (0.093  inches)  in 
diameter  and  incorporated  Kulite  XB-093-25  transducers.  As 
reported  in  Ref.  6,  the  outputs  of  both  probes  as  a  function 
of  yaw  angle  could  be  represented  well  by  polynomials  of  fourth 
order.  It  was  also  found  during  calibration  that  the  yaw 
angle  at  which  the  probe  indicated  ambient  (static)  pressure 
was  independent  of  Mach  number.  This  observation  led  to  an 
analytical  reduction  of  the  calibrations  to  simple  polynomial 
equations.  Because  the  symmetry  of  the  pressure  characteristic 
was  to  be  used  in  the  data  reduction  it  was  concluded  at  this 
stage  that  an  accurately  formed  or  machined  tip  geometry  was 
required  for  both  probes.  By  recasting  the  raw  calibration 
data  into  the  required  form,  one  complete  set  of  test  data  from 
the  compressor  was  reduced  to  velocity  profiles  and  is  pre¬ 
sented  herein.  Failure  of  the  Type  A  probe  during  handling 
prevented  further  work  with  the  first  generation  system.  The 
results  obtained  with  the  first  generation  probes  will  be  dis¬ 
cussed  qualitatively  only  since  it  was  not  possible  to  repeat 
the  early  calibration  on  which  the  reduction  was  based. 

The  second  generation  probes  were  smaller  and  the  in¬ 
clined  angle  of  the  Type  B  probe  tip  to  the  social  direction 
was  reduced  from  55°  to  35°.  The  larger  angle  had  been  chosen 
initially  so  that  the  difference  between  Type  A  and  Type  B 
probe  readings  would  be  comparable  to  the  dynamic  pressure. 

The  angle  was  reduced  because  at  55°  the  output  of  the  Type  B 
probe  became  insensitive  to  pitch  angle  at  only  a  few  degrees 
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pr-ce  tip  diameter  was  i.57  mm 


(0.062  inches)  and  used  a  Kulite  "Type  B"  screen  with  eight 
holes  as  shown  in  Fig.  8,  The  change  in  tip  geometry  (the 
larger  probes  had  a  slightly  recessed  screen)  gave  a  change 
in  the  shape  of  the  output-angle  characteristic.  It  was  also 
found  that  partial  blockage  of  the  holes  in  the  screen  (Fig. 

3)  caused  "skewing"  of  the  output-angle  characteristic,  and 
this  was  remedied  before  a  full  calibration  was  made. 

The  second  generation  probes  were  calibrated  in  a  free 
jet  and  calibration  surfaces  required  for  the  data  reduction 
using  Method  3  were  generated.  The  calibration  representa¬ 
tion  was  then  verified  by  separate  free  jet  tests  and  the 
probes  were  installed  on  the  compressor.  The  procedures  used 
with  the  second  generation  probes  will  be  described  here,  and 
the  results  obtained  in  the  compressor  will  be  analyzed  for 
their  probable  accuracy. 

Probe  Calibration  and  Verification 

A  10.16  cm  (4  inch)  diameter  continuous  free  air  jet 
exhausting  to  atmosphere  was  used  for  both  calibration  and 
verification  tests.  The  apparatus  is  described  in  Ref.  8  and 
details  of  the  procedures  are  given  in  Ref.  10.  The  probes 
were  calibrated  separately.  The  probe  mount  allowed  complete 
rotation  in  yaw  and  variation  in  pitch  from  -45°  to  +45°. 

Each  calibration  involved  recording  the  probe  output  digitally 
with  reference  data  as  the  probe  was  continuously  yawed  from 
-80°  to  +80°  and  back  at  each  of  9  discrete  pitch  angles  (-15° 
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to  +25°  in  increments  of  5°)  at  each  of  6  Mach  numbers  (0.2  to 
0.7  in  increments  of  0.1).  For  each  of  the  54  combinations  of 


pitch  angle  and  Mach  number,  300  values  were  obtained  to  de¬ 
scribe  the  probe  output  (p)  vs  yaw  angle  (a)  characteristic  for 
each  probe.  Figure  9  shows  examples  of  the  calibration  data 
obtained  with  the  Type  A  and  Type  B  probes.  The  average  of  the 
sweeps  to  and  from  was  taken  after  this  was  verified  to  give  a 
correct  result. 


In  order  to  represent  the  calibration  data  in  the  form 
required  by  Method  3,  the  following  procedure  was  followed  for 
each  value  of  Mach  number  and  pitch  angle: 


(i)  the  maximum  value  of  the  Type  A  probe  output,  pA  max» 
was  found  using  a  4th  order  polynomial  curve  fit  to 
PA  vs  aA  for  -20°  <  at  <  20° 

(ii)  the  value  of  pA  at  a* -63°  (on  the  left  branch  and 
designated  PSAL)  was  found  by  second  order  interpola¬ 
tion  over  eight  data  points 


(iii) 


the  corresponding  value  PSAR  on  the  right  branch  at 
a  =*  +63°  was  found  using  the  same  procedure 


(iv)  the  average  of  PSAL  and  PSAR  was  calculated;  pgA 

(v)  the  maximum  value  of  the  Type  B  probe  output,  pB  may 
was  found  using  a  4th  order  polynomial  curve  fit  to 
pQ  vs  a  for  -30°  <  a  <  +30° 

(vi)  the  3  pressure  coefficients,  5,  Y  and  5,  were  calcu¬ 
lated  using  the  definitions 


3  *  max  ~  PSA 
^A  max 

Y  m  max  ~  Pb  max 
^A  max  ”  ^SA 


(1) 

(2) 


and  6  ■  0  •  y 


(3) 
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(vii)  the  dimensionless  velocity  X,  defined  as  the  ratio  of 
velocity  V,  to  "limiting"  or  stagnation  velocity  ,Vt, 
was  calculated  from  the  relationship 


where 


and 


X2 


m2 


1  + 


M 


/ZcTT 


P  t 


(4) 


(5) 

(6) 


where  c^  is  the  specific  heat  at  constant  pressure 
and  Tt  is  stagnation  temperature. 

From  these  data,  surface  approximations  of  the  calibra¬ 
tion  were  obtained  in  the  form 


X  -  X ( 8, y) 


<P  *  *(6,y) 


6(j-l)' 


,  (i-1) 


(i-1) 


(7) 

(8) 


using  subroutines  given  in  Ref.  9.  The  orders  of  the  ap¬ 
proximation  were  selected  as  those  giving  the  least  error  in 
approximating  the  data  in  the  desired  range.  The  results  for 
the  second  generation  probes  are  shown  in  Fig.  10.  In  effect, 
the  calibration  data  are  completely  replaced  by  the  coeffi¬ 
cients  In  reducing  data  from  the  compressor  (or 

verification  tests) ,  the  coefficients  are  entered  into  the 
reduction  program  and  the  values  of  X  and  $  are  obtained  from 
measured  values  of  8  and  y  using  Eq.  (7)  and  Eq.  (8) . 
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The  accuracy  of  the  reduction  method  and  surface  approxi¬ 
mation  was  calculated  for  the  data  obtained  in  the  calibration 
tests.  The  results  are  given  in  Table  I  for  the  range  of  in¬ 
terest.  The  errors  are  defined  as 


ex  = 


X  -  x„ 

m  c 


100 


m 


and 


e,  =  <j>  -  6 

<t>  vm  yc 


(9) 


(10) 


where  subscript  m  indicates  the  measured  value  and  subscript 
c  denotes  the  value  given  by  Eq.  (7)  or  Eq.  (8)  using  the  se¬ 
lected  coefficients,  included  in  Table  I  is  the  error  in  the 
derived  yaw  angle  given  by 


a  —  a 
m  c 


(11) 


where  am  was  0  and  ac  was  the  yaw  angle  of  the  Type  B  probe 


corresponding  to  the  maximum  output,  p0 


max 


The  yaw  angle 


was  defined  for  the  compressor  data  in  the  same  way. 

To  evaluate  the  probable  error  involved  in  applying  the 
calibration  to  compressor  test  data,  specific  test  conditions 
were  established  in  the  free  jet  and  data  were  acquired  from 
the  two  probes  at  the  same  time  following  the  procedure  to  be 
used  in  the  compressor.  Data  were  taken  with  the  probes  set 
at  nine  different  yaw  angles  between  -65°  and  +65°.  The  data 
were  reduced  using  procedures  similar  to  those  described  for 
the  compressor  tests.  The  errors  between  the  known  (measured, 
subscript  m)  condition  and  that  obtained  using  the  measured 
data  in  Eq.  (7)  and  Eq.  (8)  (calculated,  subscript  c) ,  were 
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evaluated  using  the  definitions  in  Eq.  (9)  and  Eq.  (IQ).  The 
results  sure  shown  in  Table  II.  It  can  be  seen  that  in  a 
steady  uniform  flow,  the  maximum  uncertainty  in  the  magnitude 
of  the  dimensionless  velocity  was  1.5%,  in  yaw  angle  Q.8  de¬ 
grees,  and  in  pitch  angle  1.0  degrees. 

Compressor  Measurements 

Compressor :  The  transonic  compressor  in  which  the  mea¬ 
surements  were  made  is  a  single  stage  machine  with  axial  inlet 
flow  to  the  rotor.  Figure  1  shows  the  flow  path  to  scale. 

The  outside  case  wall  diameter  is  27.94  cm  (.11  inches)  and  the 
blade  height  is  approximately  6.99  cm  (2.75  inches),  4.83  cm 
(1.9  inches)  and  4.57  cm  (1.8  inches)  at  stations  upstream  of 
the  rotor,  downstream  of  the  rotor  and  downstream  of  the  sta¬ 
tor  respectively.  The  compressor  design  is  described  in  Ref. 

3  and  the  test  rig  in  Ref.  11.  The  compressor  is  driven  by  a 
cold-air  turbine  supplied  from  a  continuous  laboratory  air 
supply.  Instrumentation  is  provided  to  measure  stage  perfor¬ 
mance.  Radial  survey  probes  are  used  ahead  of  and  between 
blade  rows  and  the  case  wall  can  be  rotated  to  position  in¬ 
strumentation  peripherally.  A  combination  temperature- 
pneumatic  probe  of  the  type  shown  in  Fig.  11  is  positioned 
downstream  of  the  rotor  and  is  used  to  obtain  the  time  aver¬ 
age  velocity  vector  and  temperature  through  the  calibration 
procedure  given  in  Ref.  10.  Data  acquisition,  recording  and 
reduction  uses  the  system  shown  in  Fig.  4.  The  software  is 
described  in  Refs.  12  and  18. 
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Pneumatic  pressures  are  recorded  using  a  Scanivalve, 

strain  gauge  transducer,  signal  conditioning  and  DVM  to  an 

2 

uncertainty  of  ±25  N/m  (0.1  in.  of  water).  Kulite  differ¬ 
ential  transducer  signal  conditioning  circuits  output  through 
a  scanner  to  the  DVM  and  through  DATEL  amplifiers  to  separate 
channels  of  a  ±1.0  volt  range  HP  561Q  A/D  converter.  On-line 
calibration  ensures  an  uncertainty  in  Kulite  pressure  data 

to  ±0.002  of  the  pressure  equivalent  of  the  1  volt  range, 

2 

typically  ±50  N/m  (0.2  inches  of  water). 

The  geometry  of  the  rotor  and  stator  blading  at  the  tip 
is  shown  to  scale  in  Fig.  12.  Also  shown  to  scale  is  a  radial 
view  of  the  probe  (either  Type  A  or  Type  B)  and  the  range  of 
yaw  angles  within  which  it  is  rotated  in  the  data  acquisition 
process.  Velocity  diagrams  typical  of  the  flow  within  and 
outside  the  rotor  wake  are  shown.  It  is  noted  that  the  rotor 
blade  profile  is  flat  on  the  pressure  side  and  is  a  circular 
arc  on  the  suction  side.  The  stator  blades  are  double  circu¬ 
lar  arc. 

Test  Procedure:  The  compressor  speed  and  throttle  were 
adjusted  to  produce  the  desired  operating  condition  and  steady- 
state  data  were  recorded.  The  combination  probe.  Type  A  probe 
and  Type  B  probe,  with  tips  at  the  same  radial  position,  were 
adjusted  in  yaw  to  the  time  average  yaw  angle  given  by  the 
combination  probe.  Data  were  recorded  for  the  purpose  of 
calibrating  the  Kulite  probes  on-line  since  they  were  known 
to  be  sensitive  to  temperature.  This  required  recording  the 
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Kulite  probe  outputs  using  an  integrating  voltmeter  when  four 
different  levels  of  air  pressure  were  applied  to  the  reference 
side  of  the  transducers.  The  slopes  of  the  transducers  were 
calculated  using  a  least  squares  linear  fit  to  the  five  data 
points.  This  procedure  was  repeated  later,  after  the  data 
were  acquired.  The  data  were  accepted  only  if  the  two  values 
of  slope  were  within  ±1.5%  and  mean  values  were  used  in  the 
reduction.  Paced  data  at  256  positions  (counts)  across  two 
blade  spaces  were  recorded  using  a  programmed  number  (10  to 
40)  of  samples  at  each  position.  Reference  data  for  the  com¬ 
pressor  and  the  time-averaging  combination  probe,  and  the  time 
average  output  of  the  Type  A  and  Type  B  probes  were  recorded 
from  the  integrating  DVM.  The  Paced  data  and  reference  data 
recording  were  repeated  with  the  Type  A  and  Type  B  probes 
yawed  in  unison  to  four  angles  in  each  direction.  Angles 
were  chosen  to  range  from  -70°  to  +80°  with  respect  to  the 
time  average  angle  to  accommodate  the  expected  flow  yaw  angle 
behavior  within  the  rotor  wakes.  Further  details  are  given 
in  Ref .  8 . 

Data  Reduction:  The  slopes  of  the  transducers  were 
established  during  the  data  taking  procedure.  The  intercepts 
were  obtained  for  the  two  probes  using  the  time-average  mea¬ 
surements.  First,  from  the  calibration  data  for  the  probes 
such  as  are  shown  in  Fig.  9,  polynomial  (surface)  expressions 
had  been  derived  for 


Paq  *  ps 
Pt  ~  ps 


(12) 
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'.which  was 


*-p>  vi  r*  +  t  *  L\  *-*-» 


4  -3  -' - ,  —  ^  * 


bly  on  Mach  number)  and 


,  Pbo  "  ps 

-pg  =  cOg  (X,^)  =  ■  'Jj—  1  r~ 
^B0  iS0  pt  p= 


(13) 


Here,  Cp^  and  cpBo  are  pressure  coefficients  at  zero  yaw  de¬ 
fined  for  the  Type  A  and  Type  B  probes  respectively  in  terms 
of  the  corresponding  probe  pressures  pAg  and  pBg,  stagnation 
pressure  pfc  and  static  pressure,  ps-  Using  the  time-averaged 
values  for  $,  X,  pfc  and  ps  obtained  from  the  combination  probe 
in  Eg.  (12)  and  Eq.  (13) ,  values  of  pAg  and  pB^  were  calcu¬ 
lated.  It  was  then  assumed  that  these  values  corresponded  to 
the  maxima  in  the  time-averaged  output  voltages  obtained  as  a 
function  of  probe  yaw  angle.  Thus  the  two  intercepts  were 
catermined,  and  all  paced-data  samples  from  both  probes  could 
be  reduced  to  absolute  pressures. 

Data  obtained  at  each  b lade- to-b lade  location  with  the 
probes  at  nine  different  yaw  angles  were  reduced  as  a  set. 

The  same  procedure  was  carried  out  for  each  of  256  sets. 

First,  fourth  order  curves  were  fitted  to  the  curves  of  pres¬ 
sure  vs  probe  angle  (Fig.  7) .  For  the  Type  A  probe,  the  value 
of  pressure  pg^  was  found  at  which  126°  separated  left  and 
right  branches.  The  maximum  value,  PAmax»  was  obtained  at  the 
mid-point  between  the  left  and  right  branches.  From  the  curve 
for  the  Type  B  probe,  the  maximum  pressure,  p_  ,  was  de- 
rived  from  the  curve-fit  and  the  flow  yaw  angle  was  obtained 
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at  the  maximum.  The  values  of  8  and  y  were  obtained  using 
Eq.  (1)  and  Eq.  (2)  ,  and  X  and  <j>  were  calculated  using  the 
calibration  coefficients  in  Eq.  (7)  and  Eq.  (8) .  Mach  number 
was  obtained  from  a  solution  of  Eq.  (4) . 

Results  and  Discussion 

First  Generation  Probe  System;  Results  were  obtained 
using  the  first  generation  probe  system  with  the  compressor 
operating  at  15,230  RPM  (50%  of  design  speed)  and  throttled 
to  near  peak  efficiency.  A  partial  map  for  the  stage  is 
shown  in  Fig.  13,  and  the  machine  condition  for  the  first 
probe  data  is  identified  there  as  RUN  117.  Data  were  ob¬ 
tained  at  eight  radial  positions  from  hub  to  tip  across  two 
blade  passages.  The  resulting  distributions  of  dimerui  .'.nless 
velocity  X  (defined  in  Eq.  5) ,  yaw  angle,  a  and  pitch  angle, 

4>  are  shown  in  Fig.  14.  (It  is  noted  that  the  limiting 
velocity  at  room  temperature  is  approximately  760  m/sec 
(2500  ft/sec).)  Also  shown  in  each  distribution  is  the 
time-averaged  value  obtained  from  the  combination  probe. 

In  discussing  the  results,  it  is  first  noted  that  the 
data  are  shown  plotted  by  the  computer  (pen  down)  without 
smoothing.  From  oscilloscope  observations  of  probe  outputs, 
the  area  of  the  rotor  blade  wake  (which  corresponds  closely 
to  where  the  yaw  angle  is  increased)  was  clearly  unsteady. 

In  acquiring  the  data,  40  separate  samples  were  taken  at 
each  point  through  the  wakes  and  15  were  taken  outside  the 
wakes.  The  resulting  ensemble  averages  are  seen  to  result 
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in  a  smooth  and  coherent  distribution  of  flow  characteristics 
throughout.  It  is  noted  that  each  data  point  shown  was  re¬ 
duced  from  two  sets  (2  probes)  of  9  ensemble  averages  (at 
each  probe  angle)  digitized  at  quite  different  times  in  the 
acquisition  process.  The  coherence  of  the  behavior  from 
point  to  point  at  a  single  radial  position  is  therefore  bet¬ 
ter  than  might  have  been  expected.  The  clearly  coherent 
changes  in  behavior  observed  with  radial  displacement  is 
striking. 

The  agreement  of  the  measurements  with  the  time-average 
values  is  to  be  expected  since  the  on-line  calibration  pro¬ 
cedure  forces  the  time  averaged  behavior  of  the  Kulite  system 
to  duplicate  that  of  the  combination  probe.  However,  it  is 
noted  that  this  simply  fixes  the  intercepts  for  the  trans¬ 
ducers  and  the  paced  data  from  which  the  data  in  Fig.  14  was 
reduced  was  otherwise  obtained  independently. 

Qualitatively,  the  yaw  angle  distribution  is  as  would 
be  expected,  with  increased  magnitudes  resulting  from  lower 
relative  velocity  in  blade  wake  regions.  The  absolute  velo¬ 
city  magnitude  appears  to  vary  outside  the  wake  nearer  to  the 
hub,  with  larger  magnitudes  from  the  suction  side  of  the 
blades.  The  outermost  distribution  was  taken  4.83  mm  (0.19 
inches)  from  the  case  wall  and  was  probably  influenced  by 
case  wall  boundary  layer  and  tip  flow  effects.  This  is  most 
evident  in  the  derived  pitch  angle  which  is  seen  to  vary  as 
much  outside  as  inside  the  blade  wakes. 


26 


The  qualitative  behavior  of  the  results  was  examined 
using  the  first  successful  data  obtained  on  the  annulus  cen¬ 
terline  in  an  earlier  test.  The  machine  condition  is  indi¬ 
cated  as  RUN  116  in  Fig.  13.  The  dimensionless  velocity  and 
yaw  angle  results  are  shown  in  Fig.  15  and  three  discrete 
conditions  to  be  examined  at  counts  47,  88  and  95  respectively 
are  indicated.  Count  47  is  outside  the  wake,  count  88  cor¬ 
responds  to  the  maximum  yaw  angle  Cwhich  occurs  at  a  velocity 
magnitude  greater  than  the  minimum) .  The  velocity  diagrams 
which  correspond  to  the  data  at  these  three  positions  are 
shown  in  Fig.  16.  It  can  be  seen  that  the  relative  flow 
velocity  vectors  for  the  three  positions  and  for  the  time 
averaged  condition  are  at  51°  or  52°  (50.8°  was  used  in  design). 
Hence  a  relative  velocity  which  is  almost  constant  in  angle  but 
wake- like  in  magnitude  is  consistent  with  the  reduced  data.  It 
can  also  be  seen  in  Fig.  16  why  the  minimum  in  absolute  velocity 
does  not  occur  at  the  maximum  yaw  angle.  If  the  relative  flow 
angle  were  truly  constant,  then  the  absolute  velocity  vector 
would  move  along  the  broken  line  drawn  parallel  to  this  di¬ 
rection.  The  minimum  velocity  magnitude  occurs  where  the  ab¬ 
solute  velocity  is  normal  to  the  broken  line  (approximately 
true  here  at  count  88),  whereas  the  maximum  yaw  angle  occurs 
where  the  magnitude  of  the  relative  flow  velocity  is  a  mini¬ 
mum.  Hence  the  behavior  of  the  derived  velocity  magnitude 
and  angle  is  consistent  with  a  nearly  constant  relative 
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flow  angle.  The  calculated  distribution  of  the  relative  flow 
vector  is  shown  in  Fig.  17  for  the  data  from  RUN  116. 

Second  Generation  Probe  System:  Data  were  obtained 
using  the  second  generation  probe  system  on  the  centerline  of 
the  compressor  annulus,  at  18300  RPM  (60%  of  design  speed)  at 
different  compressor  flow  rates.  The  machine  run  conditions 
were  those  indicated  in  Fig.  13.  RUN  123  was  taken  near  peak 
efficiency  while  RUN  125  was  at  more  throttled  conditions. 
Ensemble  averages  were  taken  over  40  samples  at  all  blade- to- 
blade  positions.  The  results  are  shown  in  Fig.  18. 

It  can  be  seen  that  as  the  compressor  was  throttled  the 
average  yaw  angle  increased,  the  change  in  yaw  angle  through 
the  wake  increased,  and  the  wake  width  increased.  There  was 
also  a  significant  change  in  pitch  angle  in  the  wake,  with  up 
to  12°  positive  increase  indicated  at  the  throttled  condition. 
The  pitch  angle  distribution  outside  the  wake  appeared  to  be 
little  changed  with  throttling  but  showed  greater  variations 
than  was  observed  at  50%  speed.  The  velocity  magnitude  showed 
considerable  change  in  the  wake  region;  however,  examination 
also  showed  that  the  variations  in  yaw  angle  and  velocity 
magnitude  remained  consistent  with  almost  constant  relative 
flow  angle  for  each  set  of  data. 

Verification  of  the  Measurements:  The  qualitative  be¬ 
havior  of  the  data  has  been  discussed  and  has  left  few  obvious 
questions.  The  question  of  quantitative  accuracy  remains. 
Fortunately,  inherent  in  the  measurement  technique  is 
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measurement  redundancy,  which  allows  the  accuracy  to  be  veri¬ 
fied.  The  verification  at  a  given  blade-to-blade  position 
involves  an  examination  of  whether  the  data  obtained  there  at 
the  9  yaw  ancles  follow  the  same  characteristic  obtained  during 
calibration  in  steady  flow  at  the  corresponding  Mach  number 
and  pitch  angle.  Such  a  comparison  would  be  quite  simple  if 
Method  1  were  used,  since  the  pressure  characteristics  are 
there  represented  analytically.  Using  Method  3,  with  no  such 
analytical  representation  it  is  more  difficult  and  has  been 
done  properly  only  for  the  Type  A  probe. 

If  the  output  of  the  Type  A  probe  is  represented  as  a 
pressure  coefficient,  cpA,  defined  as 


it  is  found  that  cpA  depends  very  little  on  Mach  number  gen¬ 
erally,  and  very  little  on  pitch  angle  in  the  range  -5°  <  <t>  <  20° 
Thus,  when  the  local  values  of  X  (or  Mach  number)  and  $  have 
been  derived,  the  local  values  of  pfc  and  pg  are  given  by  Eq. 

(12)  and  Eq.  (13)  and  values  for  cpA  can  be  calculated  using 
Eq.  (14)  for  each  probe  yaw  angle.  These  values  can  be  im¬ 
posed  on  the  curve  of  cPa  vs  a  obtained  at  a  similar  Mach 
number  and  similar  pitch  angle  during  calibration. 

This  has  been  done  for  data  obtained  in  RUM  123  and  the 
results  are  shown  in  Fig.  19.  First,  it  can  be  seen  that  the 
time-averaged  data  obtained  as  a  function  of  yaw  angle  closely 
approximated  the  behavior  in  steady  flow.  Also,  data  acquired 
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outside  the  wake  region  (at  count  128)  was  also  consistent 
with  steady-flow  behavior.  Inside  the  wake,  however  it  is 
clear  that  the  behavior  of  the  ensemble-averaged  data  do  not 
follow  the  characteristic  established  in  steady  flow  except  at 
the  very  center  (at  count  188).  Going  into  the  wake  (at  count 
181) ,  the  characteristic  is  skewed  and  the  peak  probe  output 
is  to  the  right.  Coming  out  of  the  wake  (at  count  195) ,  the 
characteristic  is  skewed,  with  the  peak  output  to  the  left. 

It  must  therefore  be  concluded  that  reduced  data  ob¬ 
tained  outside  the  unsteady  wake  regions  can  be  accepted  to 
within  the  uncertainty  derived  from  the  measurement  uncer¬ 
tainty.  The  reduced  data  for  the  area  of  the  wake  contains 
an  unknown  uncertainty  for  which  an  explanation  is  required. 
The  values  at  the  center  of  the  wake  however,  appear  also  to 
be  verified. 

Assumptions  and  Sources  of  Error:  The  verification 
suggests  that  for  the  measurements  made  outside  the  wake, 
the  method  gives  acceptable  results.  Several  assumptions 
were  made  nevertheless,  which  should  be  examined. 

The  assumption  that  the  probes  respond  with  negligible 
error  to  fluctuations  at  blade  passing  frequencies  was  sup¬ 
ported  by  oscilloscope  observations  of  raw  signals.  The 
probes  were  observed  to  respond  quite  differently  in  the 
wake  region  compared  to  outside  the  wake  region.  The  probe 
output  fluctuated  only  over  a  clear  fraction  of  each  blade 
passage  and  was  otherwise  quite  steady  when  triggered  by 
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the  one-per-re”-  signal.  That  the  flow  outside  the  wake 
was,  on  the  average,  steady  in  the  rotor  frame  was  clearly 
demonstrated  by  the  same  observations. 

The  assumption  was  made  for  the  on-line  calibration 
that  the  time  averaged  behavior  of  the  Kulite  probe  system 
could  be  equated  to  that  given  by  the  combination  probe. 

This  was  supported  in  the  verification  check  above,  however 
the  result  was  that  the  absolute  levels  of  the  parameters 
derived  from  the  Kulite  probe  system  depended  on  the  ac¬ 
curacy  of  the  combination  probe.  Since  averaging  errors 
can  occur  in  pneumatic  probes  (although  as  shown  in  Ref.  10 
they  are  negligible  for  the  conditions  found  behind  the 
rotor  at  present  speeds) ,  errors  in  the  pressure  measure¬ 
ments  of  the  combination  probe  sensors  were  assumed  and  the 
effect  on  the  Kulite  measurements  was  evaluated.  It  was 
found  that  assumed  errors  of  as  much  as  10%  in  the  pneumatic 
pressures  gave  rise  to  negligible  differences  in  the  velocity 
and  yaw  angle  results,  and  differences  ranging  from  -h°  to 
+1°  in  pitch  angle  results. 

The  source  of  the  discrepancy  in  the  behavior  of  the 
probes  in  the  wake  regions  is  not  resolved.  It  is  encouraging 
that  the  measurements  near  the  wake  center  may  be  acceptable 
since  the  knowledge  of  the  wake  width,  the  minimum  velocity 
and  the  conditions  at  the  wake  edge  are  sufficient  for  many 
engineering  purposes.  An  attempt  will  be  made  however  to 
resolve  the  remaining  uncertainty  in  the  method.  Of  the 
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possible  explanations,  the  suggestion  that  there  is  a  vary¬ 
ing  temperature  effect  through  the  wake  has  been  rejected 
since  the  thermal  lag  of  the  transducer  sensor  could  not 
allow  such  a  rapid  response.  The  remaining  potential  expla¬ 
nations  include: 

(i)  the  steady  flow  error  of  probe  measurement  in  a 
shear  layer 

(ii)  differences  in  unsteady  response  at  different  yaw 
settings 

(iii)  incorrect  averaging  procedure  in  a  flow  in  which 
Mach  number  and  flow  angle  are  unsteady  and  probe 
system  response  is  non-linearly  dependent  on  both. 

A  limited  attempt  has  been  made  to  investigate  errors  which 
might  stem  from  the  averaging  procedure.  However  it  is  dif¬ 
ficult  to  see  why  a  lack  of  symmetry  in  the  probe  output 
could  be  caused  this  way.  It  is  noted  that  the  error  which 
would  normally  result  from  probe  measurement  in  a  steady 
shear  layer  is  consistent  with  the  behavior  observed  quanti¬ 
tatively,  and  will  be  investigated  further.  The  possibility 
of  differences  in  unsteady  response  at  different  yaw  angles 
is  to  be  investigated  in  the  near  future  by  changing  probe 
tip  geometry. 

IV.  CONCLUSIONS 

The  Dual  Probe  Digital  Sampling  (DPDS)  was  applied  to 
measure  the  blade-to-blade  and  hub-to-tip  flow  field  from 
two  blade  passages  of  a  compressor  rotor.  All  three  com¬ 
ponents  of  the  velocity  were  derived  with  relative  success. 
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The  resuica  showed  a  coherent  and  realistic  behavior  near 

peak  efficiency  at  5Q%  speed.  Broadening  of  the  wake  with 
increased  yaw  and  pitch  angles  was  measured  at  mid  span  and 
the  compressor  was  throttled  at  60%  speed.  Outside  the  wake 
regions,  the  flow  was  found  to  be  steady  in  the  rotor  frame 
and  the  data  were  shown  through  inherent  measurement  redun¬ 
dancy  to  be  valid.  Inside  the  wake  regions,  the  flow  was 
found  to  be  unsteady  in  the  rotor  frame  and  the  reduced 
data  contain  an  unknown  uncertainty  except  near  the  wake 
center  where  they  are  thought  to  be  valid.  Attempts  are  to 
be  made  to  resolve  the  remaining  uncertainty  and  to  compare 
the  measurements  with  those  from  an  LDA  system. 


33 


TABLE  I.  Errors  in  Surface  Approximation 
of  Probe  Calibration  Data 


Test 

Condition 

Approximation 

Error 

X 

$  (deg) 

ex 

(deg) 

ea  (deg) 

1347 

15.00 

-1.3357 

.17 

-  .69 

1349 

10.00 

-  .6396 

.04 

-  .52 

1346 

5.00 

.2487 

-  .03 

-  .15 

1345 

0.00 

1.0161 

-  .07 

-  .40 

1750 

15.00 

-  .9074 

.25 

-  .73 

.1758 

10.00 

.1439 

.13 

-  .71 

,1750 

5.00 

.  .1923 

-  .33 

-  .68 

.1751 

0.00 

.0256 

-  .22 

-  .53 

.2163 

15.00 

.3642 

.96 

-  .13 

.2161 

10.00 

-  .0087 

.91 

-  .64 

.2160 

5.00 

-  .8027 

.97 

-  .80 

.2162 

0.00 

-  .9259 

1.36 

-1.05 

.2626 

15.00 

-  .6210 

-  .36 

-  .16 

.2624 

10.00 

-  .4118 

-  .56 

-  .08 

.2622 

5.00 

-  .6578 

.98 

-  .19 

.2622 

0.00 

-  .8000 

-  .54 

-  .91 

.2944 

15.00 

-  .5198 

-  .20 

-1.01 

.2945 

10.00 

-  .8086 

-  .64 

-  .38 

.2945 

5.00 

-1.1207 

-  .03 

-  .76 

.2948 

0.00 

.1303 

-  .41 

-1.16 
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TABLE  II.  Errors  Obtained  in  Verification 
Tests  in  Steady  Flow 


Test 

Condition 

Total 

Measurement 

Error 

X 

<t> 

(deg)  a 

(deg) 

£X  (%> 

ea  (deg 

1420 

0 

0 

-0.04 

-0.04 

-0.31 

1459 

5 

0 

+  1.49 

-0.26 

+0.11 

1459 

5 

0 

+  1.52 

+0.21 

-0.1 

1450 

10 

0 

+  1.28 

-0.18 

-0.24 

1448 

15 

0 

-0.02 

+0.18 

+  0.30 

170 

10 

0 

+  1.22 

+  0.22 

-0.78 

170 

10 

5 

-1.17 

+  1.03 

-0.45 

170 

10 

10 

-0.31 

-0.2 

-0.22 

2049 

7 

0 

+  0.53 

-0.66 

+0.18 

211 

10 

0 

-0.31 

-0.16 

-0.2 

2052 

13 

0 

+  0 . 26 

-0.21 

-0.24 

> 
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Figure  11.  Views  of  the  Combination  Probe 
With  Thermocouple  Sensor  Removed.  Right  -  Thermocouple  Sensor) 


0.00  64.00  128.00  192.00  256.00 


/ 

/  c 

/  I.  1  . 1 

DmENSIONLESC  VFL?C!TV  SCALE 


Fig.  16.  Velocity  Diagrams  for  Different 
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